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^_ Abstract: We investigate the role of the choice of the upper phase space 

limit Q in the Curci-Furmanski-Petronzio (CFP) factorization scheme, which 
exploits dimensional regularization MS scheme. We examine how the 
choice of Q influences the evaluation of the standard DGLAP (inclusive) 
evolution kernels, gaining experience needed in the construction of the ex- 
'■j~ H elusive Monte Carlo modelling of the NLO DGLAP evolution. In particular, 

O | we uncover three types of mechanisms which assure the independence on Q 

of the inclusive DGLAP kernels calculated in the CFP scheme. Wc use the 
examples of three types of the Feynman diagrams to illustrate our analysis. 
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1. Introduction 

This study is part of an effort aiming at construction of an exclusive 
Monte Carlo modeling of DGLAP pQ evolution of the parton distribution 
functions (PDFs) in the next-to-Leading-Order (NLO) approximation using 
work of Curci-Furmanski-Pertonzio (CFP) [2] as a starting point. Standard 
inclusive PDFs are defined within the framework of the collinear factoriza- 
tion theorems El [5]- The ongoing project of defining and implementing 
in the MC form exclusive PDFs (ePDFs), see ref. [6j [7], sometimes also 
referred to as fully unintegrated PDFs [8] , is based on the older formulation 
of the collinear factorization of ref. [3] reformulated later on by CFP [2]. 
The CFP work uses dimensional regularization in MS scheme and physical 
axial gauge. 

The construction of ePDFs in the Monte Carlo form requires defining 
and calculating new exclusive evolution kernels. Moreover, it is critical 
to understand and analyse the properties of exclusive kernels, especially 
cancellations of the infrared singularities diagram by diagram due to gauge 
invariance, see study in ref. [S]. 

In this contribution we comment on the issue of the independence of 
inclusive/exclusive kernels on the choice of the upper phase space limit Q 
in their evaluation based on the Feynman diagrams. Of course, the inde- 
pendence of the inclusive DGLAP evolution kernels in the CFP may be 
regarded as obvious and trivial. However, in the actual calculation of the 
kernel from the Feynman diagrams the genuine mechanism which protects 
its independence on Q looks rather mysterious and not obvious at all. The 
choice of Q turns out to be important in the analytical evaluation of the 
NLO kernels in the CFP scheme, because it determines quite rigidly the 
parametrization of the two-parton phase space. In addition, in the con- 
struction of Monte Carlo model for ePDF the same upper phase space limit 
variable Q is closely related to the evolution time variable. It is therefore 
quite interesting to have a closer look into the above phenomena. 

In the following we will show that there are three different mechanisms 
which assure the independence of inclusive NLO DGLAP evolution kernels 
on the upper phase space limit Q in the CFP scheme. We demonstrate each 
mechanism using an example of the Feynman diagram contributing to NLO 
DGLAP kernel. We shall use subset of diagrams shown in Figure [T] 

2. Notation 

We consider two-gluon real emission diagrams. For the four- momentum 
parametrization we use Sudakov variables: 



h = aip + (3in + k i± , i = l,2, 



(1) 
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Fig. 1. Example Feynman diagrams contributing to NLO DGLAP kernel. 



with p being the four-momentum of the incoming quark and n a light- 
cone vector. Four- vectors of two emitted gluons are k\ and k 2 , with their 
transverse parts being ki± and k2± respectively, and k 2 = (ki + k 2 ) 2 being 
their effective mass. Since the emitted gluons are on mass shell and we are 

k 2 k 2 

in the massless theory, f3i are fixed and equal to /3j = 2«/(pn) = 2a/(pn) ■ 

We 

will also use q symbol for the off-shell momentum q = p — k = p — k\ — k 2 . 

In the CFP scheme [2] the contribution of each Feynman diagram to the 
DGLAP kernel is extracted from the phase space integral: 

P(x) = Res (Jd$6(l-x-ai- a 2 ) p(k u k 2 ) @(s(ki,k 2 ) < Q)\ , (2) 

where Reso is the residue at e = (the coefficient in front of - pole in 
the dimensional regularization) , 5(1 — x — a\ — a 2 ) is the definition of the 
Bjorken variable, p is a contribution from a Feynman diagram (originating 
from 7-traces, etc.). The element of the two gluon phase space <i$ is given 
by: 

d n h n ^d-k 22w6+{kp 



: 2n5 + (k() 



(3) 



(2vr)™ "" v " i/ (2vr) 

The theta function in equation ^ encloses the phase space from above 
using dedicated kinematical variable s(ki,k 2 ). The choice of phase space 
enclosing, s(k\,k 2 ), determines the choice of evolution time variable in the 
construction of the MC implementation of ePDF. There are many possi- 
ble choices for s(k\,k 2 ) function, here we will concentrate on two of them: 
s(ki,k 2 ) = max{|kij_|, |k2j_|}, which corresponds to the transverse momen- 
tum evolution time and s(k\,k 2 ) = max{|ai|, | a.2 1 } , which corresponds to 
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rapidity related evolution time. Scalar quantity a, is defined as a modulus 
of the vector variable: 

a, = ^ (4) 

and we call it angular scale variable. It is related to rapidity via y% = ln|aj|. 
These two cases will be respectively referred to as phase space with k±- 
ordering and angular-ordering (a-ordering) . Other popular choices of s- 
function include total virtuality \/—q 2 and maximum /c-minus, max(/cj~ , k^ )• 
Since we will show calculations in both angular ordered and fcj_-ordered 
phase space we give the phase space parametrization in both sets of variables 
(remembering that we work in dimensional regularization with number of 
dimensions equal to n = 4 + 2e, e > 0): 



d*k. = ^■^^^^^n.dh^J&l^, (5) 

and 



1 ^l+2e datl da2 2 +2e 2+2t . n A J l+2e l+2e , a \ 

d ° = 4^ (27T)6+^ V V 1 " 2 dn ^ da ^ a t a 2 ■ ( 6 ) 

3. General structure of kernel contribution 

Having in mind that we want to investigate the mechanism of the inde- 
pendence of evolution kernels on the choice of the variable s used to enclose 
phase space, let us look more closely into the phase space integral of kernel 
contribution using k± -ordering and a-ordering. 

General structure of kernel contribution is given by equation Q. The 
presence of the residue Reso ensures that only part proportional to single 
pole ~ contributes - this has to be kept in mind. For fc^-ordering the 
distribution p(ki,k 2 ) has general form: 

p(h,k 2 ) = Cg 4 1 T{k x M\ (7) 

r{ki,k 2 ) 

where C is the color factor specific for each diagram, g is related to strong 
coupling by g 2 = 2(2tt)o>s, T(k\,k 2 ) is dimensionless function and q 2 = 
^pk^ + ^^k| ± + 2kij_k 2 _i_. For a-ordering we have: 

p(ai,a 2 ) = Cg A „ 2 ^ T , (a 1 ,a 2 ), (8) 

a{a l 2 q 4 (ki,k 2 ) 

where T(ai,a 2 ) is dimensionless and q 2 = ~ a2 a| + a| + 2aia 2 . The 
above specific form of p enables immediate factorization of one e pole due to 
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integration over the overall scale variable Q, which we explicitly introduce 
by means of the identity 1 = f® dQ8(Q = s(ki, k 2 ))- The remaining integral 
is parametrized using dimensionless variables y\ = ^W* or ?/j = ^: 

(*) = Reso ( (ffi 2 C ? ^ / ^ ^ 
v 7 1 V27T/ (2vr) 2 + 4e y ai a 2 

x 5(1 - x - cti - a 2 ) [ dn^^r f dQQ u ~ l ( 9 ) 
J L 1 e Jo 

x f dy' x dy'My'2) l+2e ^444 Sil-m^y',})}, 

JO 1 ^1^2:") J 

and 

P (x) = Resoj {-) Cj-^ J — — (a^) 

X <5(1 - x - ai - a 2 ) / dn^ e 4r / Q dQQ Ae ~ l (10) 
x /' ^i^ 2 (yii/ 2 ) 1+2e T ^ l,V2,9 l S(l - max{ yi ,y 2 }) \. 

Jo (r{yi,y2,0) 



Now in eqs. feb and ( 10 ) the pole - gets explicitly factorized off in form of 
the integral f dQQ 46 ' 1 = and it is now transparent that the integrals 
of the above equations feature at least single - pole. 

Possible additional - poles may arise from internal singularities of the 
integrands of Feynman diagrams. They are always connected with integra- 
tions over transverse degrees of freedom (yi). The longitudinal components 
can also lead to infra-red (IR) singularities, when a» — > but this type of 
singularities do not lead to additional e-poles, because in the CFP scheme 
they are regularized in a non-dimensional manner]^} 

Furthermore, equations ^ and ( |10| ) show explicitly the differences be- 
tween exclusive kernel contributions (integrands). This means that exclusive 
MS evolution kernels do depend on evolution time variable. 

Equations ^ and (10) are the starting point for the investigation of 
the dependence of inclusive evolution kernels on the upper phase space 
limit/evolution time variable. There will be at least two cases to be con- 
sidered: (i) with no additional internal singularities present, hence terms 



For regularization of IR singularities CFP use principal value prescription: — — > 
a 2 < l S 2 ■ We also use the following notation of CFP for divergent integrals: 

fa da ^fs? = Io and Jo da ln a ^Ts? = 7l 
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originating form the expansion of (aia 2 ) 2e = 1 + 2eln(aia2) can be ne- 
glected, (ii) with the additional e-poles due to internal singularities present, 
hence the expansion term 2eln(aia2) is contributing. 



4. Kernels independence on the evolution time variable 

In this section we shall comment on three mechanisms which in the CFP 
factorization scheme actually protect the independence of the inclusive NLO 
DGLAP kernels of the way the phase space is closed from the above. It will 
be demonstrated using example Feynman diagrams. 



4-1- Case 1 - no internal singularities 

Here, the independence will be demonstared using the example interfer- 
ence diagram of figure [2j Starting with the expression of equation ( 10 ) the 




Fig. 2. "Crossed ladder graph", free from any internal singularities, 
calculations can be carried out in 4 dimensions 

P% x (x) = N / ——5(1 - x - ai - a 2 ) 

J ai a 2 Jo 



2tt 



x / dy\dy 2 yiy 2 



ai a 2 

T(yi,y 2 , 



(11) 



q 4 (yi,y 2 ,(f>) 



5(1- max{yi,y 2 }) , 



where is normalization constant. In massless QCD the integrand has a 
nice property that y\ and y 2 integrals can be combined together into one 
integral over the whole space by means of a simple change of variables 
VI =y,V2 = l/y, then: 



~ , > Ar fdaida 2r „ N f 2n ,, f°° , T(y,l, 

P% x (z )=N / -± 5 (l-x-a x -a 2 ) dcP dy y^f^- 
J ai! a 2 Jo Jo TKVA, 



(12) 



Since the change of the phase space enclosure from angular scale to trans- 
verse momentum translates into linear change of variables, ki± = aiai, 



2 It results from the fact that f or T depend on the ratio yi/3/2 only, i.e. T(yi, y2,4>) = 
T(\yi,\y2,4>). 
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y[ = ttiUi, and now the y integral extends from zero to oo, hence the joint 
integral is manifestly the same for both kinds of phase space enclosure. 

The argument presented in the above example holds for all kernel con- 
tributions free from internal singularities. 



4-2. Case 2 - diagram with internal singularity minus counterterm 

Second case is represented by the double bremsstrahlung diagram, see 
figure [3j It has an internal singularity when one of the emitted gluons is 
collinear (the other one kept non-collinear). The additional contribution to 
the residue due to the additional (double collinear singularity) terms is of 
a type 2eln(«ia2) x ^. This diagram is special because it is accompanied 



jvwv/ywvi 

WW\AAAA/l 



UaaaAaaaaJ 



kAAAAAAAA/V 



Fig. 3. Bremsstrahlung graph accompanied by its soft counterterm, both featuring 
double e-poles. 

by the soft counterterm, which is simply a square of leading-order (LO) di- 
agram. The soft counterterm is present due to the factorization scheme (by 
construction), see ref. [UEIITO]. We will show, that in this case, the indepen- 
dence of inclusive kernel contribution on the upper phase space limit / future 
evolution time variable is assured by the presence of the counterterm, which 
will cancel the additional ~ \n{a.\a.2) term. 

Since we have shown in section |4.1| that only terms leading to additional 
e poles can lead to differences between the two choices of evolution time 
variable now we will concentrate only on them. The singular contribution 
of the double bremsstrahlung diagram is of a form: 

pasmg^ _ j^gJ ^ f ^±^1^1 _ % _ a% _ a2 ) ( aiQ , 2 ) 2e f df^+L 

[ e J Oil a.2 J 
x ^ diftcftttto)^ - ^-—^ ^6(1 - max{ yi , y2 }) |. 

(13) 

Now combining of the two phase space integrals is not possible any mord2 



3 In n — 4 gluing two y;-mtegrals is still possible using the cut-off regularization. 
However, one has to watch out for the cut-off dependent integration's limits. 
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due to the presence of the term (yi£/2) 2<E (coming from phase space) and 
regularizing 1/yf singularity. There will be differences between k± and 
a parametrizations. There are two sources of this differences. The first 
one is simply the difference between the integrals in both parametrizations. 
The second is the mixing of double pole and e terms from the phase space 
factor (aia2) 2e . The difference between the two phase space enclosures 
(parametrizations) is 



p: 



a— ft j 
Br 



Resold f^ 2 /^^ 
e \2tt J J a\ 012 



5(1 



u\ — a2)2T 2 



In 



«2\ 



In 



— + 2em(aia 2 )- 

OL2 e 



(14) 



comes from the 



where coefficient f 2 = j(l + (1 - «i) 2 )(l + x 2 /(l - ai) 
product of the numerators of two LO kernels. 

For the counterterms, which are much simpler, due to their LO structure, 
the difference between integrals in k± and a space is only due to the phase 
space factor a| e : 



P c r^(x) = Res 




da\da25(l — x — a\ — a2)21n(a 2 ) 



2T 2 



ai«2 



(15) 



It is almost manifest that the integrals of eqs. ([14j) and (15) are the same, 
which means that 

(x) = P^(x). (16) 



r Br 



Summarizing, there is an exact cancellation of differences between the two 
choices of upper phase space limit among the bremsstrahlung diagrams and 
their soft counterterms. It results in the independence of the kernel contri- 
bution on the two choices of the phase space enclosure variable under con- 
sideration. On general ground, the same statement should hold for other 
choices, for example for the total virtuality. However, in case of virtuality it 
is much more difficult to show, without actually performing the integration, 
that the final result for the inclusive NLO kernel is the same as in the above 
cases of /cj_ -ordering or a-ordering. 

We want to emphasize the crucial role of MS'-like terms ((aia 2 ) 2e and 
a^) in restoring the independence of the kernel contribution on the choice 



of phase space encloser, which can be seen explicitly from equations (14) 



and (15). 



In the above discussion we were analyzing certain contributions from 
double bremsstrahlung diagrams and soft counterterms as representing the 
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difference between the cases of phase space enclosure using maximum k± 
or, alternatively, maximum angular scale a. In fact, these terms are com- 
pletely absent in case of maximum k±, which means that the maximum k± 
is effectively representing formal scale parameter \x of the dimensional reg- 
ularization MS. These terms are also nonzero for other popular choices of 
the phase space enclosure like maximum /c-minus and total virtuality. 

4-3. Case 3 - two-real-gluon internal singularity versus virtual diagrams 

The third case is a class of diagrams with an internal singularity due 
to parton pair emission from the ladder, for which the independence of the 
evolution time variable in assured by the corresponding virtual diagrams. 
The example diagram of this type is shown in Figure |4j This is a diagram 



Fig. 4. Gluon pair production diagram and the corresponding virtual diagram 
(vacuum polarization). 

with gluon pair production, where the internal singularity occurs when the 
invariant mass of the produced pair goes to zero. The additional e pole 
originating from the singularity 1/k 2 = l/(k\ +k2) 2 , together with (aia2) 2e 
factor, will lead to a familiar mixing terms in the residue. This mixing 
terms leads to the differences between real parton integrals once the different 
choices of the upper phase space enclosure are applied. Of course, in CFP 
scheme there is a mechanism which brings back the independence of the 
inclusive kernels on that. In this case the contributions of the corresponding 
divergent virtual diagram do this job. 

In this case the calculations are technically more complicated and have 
to remain beyond the scope of this contribution. In fact the independence 
was explicitly checked by means of switching from the angular scale to the 
overall virtuality as the upper phase space limiting variable as they are best 
suited for the singularity structure of these diagrams. 



We investigate the mechanism which ensures the independence of the 
NLO DGLAP evolution kernels calculated within Curci-Furmanski-Petronzio 
scheme on the choice of the upper phase space limiting variable s(k±, fo). 




5. Conclusions 
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It was shown that for different groups of Feynman diagrams there are 
three mechanisms which work in order to compensate the differences due to 
change of the type of s(k\, k?). The independence is demonstrated explicitly 
in case of transverse momentum and rapidity related variable a. (The inves- 
tigation has been carried out also for different choices like overall virtuality 
q 2 , maximum light-cone variable fc-minus, but no details are reported here.) 

We have show that the mechanisms protecting this property involves 
either soft counterterms of CFP scheme or virtual diagrams. 

In case of the MC implementation of the exclusive PDFs, see refs. [61 
El [10] , keeping track of these phenomena in the kernel calculations is useful 
for understanding what happens while switching from one version of the 
evolution time variable in the Monte Carlo to the other, more details will 
be provided in ref. [TT] . 
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